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CHAPTER  I 
INTRODUCTION 

Aii  automatic  method  for  the  determination  of  the  numerical  solution 
of  integro-differential-ili  fterence  equations  iu  presented.  The  basic 
form  of  an  integro-difforential-differonce  equation  is  assumed  to  be; 

Ly(t)  + J'(L-)  -h  j ::(\)  d\ 

with  initial  values  // (o)  , whore  u and  f arc  real  vector  functions 

of  the  independent  variable,  t-.  L and  G are  real  matrix  functions  of  t. 

The  numerical  method  1:;  an  extension  of  the  method  of  Frobenius , 
Lapidus  and  Sienfeld  (1971),  Kreysziy  (1967).  The  basic  theory  of  exis- 
tence and  uniqueness  for  this  equation  is  known,  and  some  basic  references 
are  [baker.  Miller,  Cushing,  and  Bellman  and  Cooke]. 

Background 

Mathematical  models  in  applied  science  are  used  to  represent  and 
simplify  complex  problems.  Mathematical  models  ire  deterroiiaistic  or 
stochastic.  Deterministic  models  include  linear  programs,  non-] inear 
prc'K;r..i«s,  differential  e.  mat  ions,  and  other  specialized  forms,  such 

as  'nhqro-difforenfial-d i : i.orence  equations.  C.T.II.  Baker  (19/4),  states 

''there  Ins  been  litU  or  no  invests;  its  •:>.  in  tlic  published  lilern- 
: ure  of  the  eld  Lot  »•  autonvitie  it.p!  . . ..  ■ ion  of  raetliods  for 
Vol terra  intogro-di:  ; . rential  equal  ion 

lVLvct  further  described  jotential  tea;  exi."  ..  for  boundary  value 

problems  in  intoyro-d  i t ! ■ l . -:’.Linl-diL£orci.c-'  .. .rations.  His  aa^virent.  tueme 

is  the  programming  dill  lculty  of  numerical  .f  .nation  of  initial  value 

problems;  and  LliaL  seems  to  Lx1  reason  cnou.f.:  to  avoid  Uie  more  dal.,  icult 


Ixaundary  value  problem. 


Soioe  integral  equations  and  integro-dif ferential  equations  are  easily 


j 

i 


I 

\ 


i 

« 


transformed  into  differential  equations.  However,  this  is  usually  diffi- 
cult for  integro-dif ference  equations  and  integro-dif ferantial-difference 
equations,  or  in possible. 

We  will  present  a method  for  the  automatic  numerical  into; ration  of 
initial  value  problems  in  integro-dif  forential -difference  equations,  and 
integro-dif foronce  equations.  This  method  cun  then  be  used  wild:  the  codes 
of  Doiron  (1%7) , Hunter  and  Childs  (1177),  and  Schuetz  (1976)  to  develop 
an  automatic  solution  of  boundary  value  problems.  These  codes  tire  '.used 
upon  tiae  power  series  integration  method,  as  is  the  code  presented  here. 

Hunter  reported  success,  using  die  power  series  method,  in  regression 
analysis  with  differential  equation  models.  Dr.  Bart  Childs  suggested 
that  power  series  integration  could  be  an  effective  bool  for  the  solution 
of  integro-dif £erential-dif ferer.ee  equations,  and  integro-dii:  f erence 
equations.  A primary  advantage  of  tec;  power  series  moth.d  is  that  it 
provides  a convenient  method  for  inclusion  of  the  "history"  or  "memory" 
in  calculations,  see  Nachlinger  aid.  Wheeler,  1973. 

Integro-dif ferential-difference  equations  arise  in  a number  of 
applications.  The  most  frequent  cure  related  to  Problems  oi.  Heat 
Conduction  - Nachlinger  and  Wheeler  (1973) , Nuclear  Transfer  and 
Transport  - Delves  and  Walsh  (1974),  Population  Models  - Cushing  (1976), 
Renewal  Theory  - Feller,  (1941),  and  V : scoelasticity  - Malone  (1971). 


t 


CHAPTER  II 


A POWER  SERIES  REPRESENI'ATIO: . OF  AN 
I NTEGR>-DIFFERf2'3CE  EQUATION 


A power  series  is  taken  tu  be  an  infinite  series  of  the  form; 

| eH>  (t-af  . c,™  . cM(t -a)  . * ... 

'i—O 

where  e ...  are  constants  called  coefficients  of  the  series,  a xs 

a constant  called  the  center  of  expansion,  and  i is  U:a  rnuependent  v,. tri- 
able. Notice:  indexes  are  parenthesized,  exponents  are  not. 

If  in  particular  a=0,  we  obuain  a poser  series. 

I cli>  t* . e<°>  (2-2 

1=0 

The  basic  idea  of  the  power  series  representation  of  integro-dif ference 
equations  is  straightforward. 

Given:  A scalar  integro-dif ference  equation: 

yd)  = J g(t-\)  yOJ  d\  (2_J 

o 

We  can  represent  the  result  of  the  integral  as  a power  series  in  t. 

fvU-X ) y(X)  d>.  - j h(i>  t*  - h<°>  . i:<l>  t . h(2>C  . ... 


where  the  h(o) , h(1) , ...  are  f nations  of  the  coefi  cients  of  the  ori  ;;  • 
series  for  j and  y.  The  limbs  of  integration  yield  = 0.  'lb  .illus- 
trate, v.e  represent  tlio  given  functions  y(\)  and  ■;  (!->■)  in  power  series: 


3 


4 


„ . Y (k)  .k 
ij(\)  = 2,  y * 

k=o 


km)  - i 9la>  u-xr 

m=o 


(2-5) 


(2-6) 


Substitutin' j these  series  into  equation  (2-4)  yields: 


j'r  y j(k)  \k)  ( y /'J  ( ->-f  d\)  = i i 

o k -D  fU—O  t—O 


(i)  i 


(2-7) 


The  results  of  the  multiplication  of  the  two  power  series  l.uy  be  rewritten 
as  expressions  giving  coefficients  of  like  powers  of  t: 


V - < 

l n 


*>t*  , l r \ y !<*-*-*>  2 a'j] 

— /n  »“  wi  -7:r_V> 


(2-3) 


•i-o  t-o  e 3~° 

Inbercluc.ging  the  order  of  independent  operators  of  surra tiun  and  integra- 
tion gives: 


i-1 


t . 


' JiKi  - l l g(i~2~j)  / A-7’  dx  (2-9) 

i=l  3=0  o 


1 n & 


'the  limits  of  integration  yield  a particularly  canpact  result: 


J'k: 


i y ! ~x>  dx  - ?wr.' 


. k+j+1 


(2-10) 


Substituting  the  results  (2-3  0)  into  equation  (2-9)  and  co  llecting  coeffi- 
cients of  lire  powers  of  t yields: 


/ -i  { ).  ')  j(t-X)  uA 
o 

■ i [ I <» 


or 


y <yh> 


. j-  ,'i-23 


'id-zy  J 


t 

(2-11) 


h(i)  = l 

t i ’ ~ ' 

We  rewrite  ox 

h(i)  - i 

ti.i '.  ice  that  tiie  li 
',\e  again  call 
superscripts  to  tic 
ci  :e  C £ ic  ients  ti  la  t 


o 


f;/j;  1~^))  [r  ‘j 


r A;  ( 

i-J  l-J+l 


■k'!. 

i-i+a-vy  J 


(2-12) 


o') 

nation  (2-12)  in  a form,  in  which  the  values  or  k are: 


v-  (J)  (i-l-j) 

) j y 


[ 


ct 

-l) 


intcgratior 


ctrve 


(o) 


(2-13) 


the  reader’  s attention  to  Lire  use  of  parenthes iteci 
mote  coefficient  indexes.  This  index  nay  apply  to 
are  scalars,  vectors,  or  tutrices. 


CHAPTER  III 


■ 


A SCALAR  .1  NTaUW-Dlt’FERtNl'lAlj-U  I KK  'J<ENCE  1 U h \T  1 1 1 -i 

When  a state  of  tlu  system  is  arbitrarily  specified  at  the*  ir..tuil 
time,  y - y(o)  at  t ••  0,  then  a solution  exiswJ  ;or  t > 0 and  is  • unruely 
determined  by  these  carnations  Richtmyer  and  ."orlon  (196/).  iiu  • • • ■ 
or.  the  scalar  integro-clif  ferential— difference  equation  is  r.bsuiw.  ‘-o 

t 

'j(t)  = f(y(t))  h J g(b-\)  y(\)  d\  ('->) 

o 

with  initial  values  y (o) . The  functions  and  j are  real  valued  ■-■-•c- 

tions  of  the  independent  variable,  t. 

Ke  denote  the  power  series  for  y as 

sin  . | /",«  <■«> 

i-o 

From  term  by  term  differentiation  of  equation  (1-2)  it  is  easily  : . lined 


y(t)  = l (i+1)  y(t+1)t°  (•' 

i=o 

Furthermore  let  f(y(t))  and  g(\)  be  represented  by  the  following  lower 


series: 


f(u(t))  = l /“V 


y(t-X)  « j g(i)  (t-X)1  (3-1) 

Wb  are  assuming  that  given  i / we  can  calculate  j and  as  we 

develop  additional  y coefficients,  we  can  develop  additional  f coefficients. 


Tha  integral  in  equ.it ion  (3-1)  my  be  repres  ntol  by  a power  • 


^ f ' ) • 

/ £?rt-x;  ya;  dx  -■  I 


(3-5) 


The  coo  i f icicnts  arc  calculated  by  the  equation  (2-12)  or  (2-13).  Again, 

' ; i 

the  limits  or  integration,  give  '■  1 ''  -0.  We  can  obtain  a :ourrence  relation 
for  tlie  higher  coefficients  in  equation  (3-1)  . By  the  proper  power  scries 
substitutions  equation  (3-1)  becomes: 


v f \ (i+1) ,i 

l (i+l)  y t 


l j 


XI)  A 

u-  -r 


r • (i)  ■ ■* 
L ,l  L 


(3-6) 


z=o 


^=o 


Coefficients  of  like  powers  of  c yields: 


(i -hi)  y(i+1)tl 


or 


y 


(l+1  ) 


/i}tl  + h^t1 


rJi>  , hH>) 


v+1 


(3-7) 


(3-S) 


(o) 


vdth  initial  values  y(o)  - y""  • Through  this  relation,  additional  coeffi- 
cients of  the  power  series  solution  can  be  obtained. 

Wo  assume  the  f(y)  in  equations  (3-1),  (3-4),  (3-7)  and  (3-3)  to  Le 

of  the  classes  of  power  series  operators  discussed  in  Appendix  A.  Tint 
(o) 


is,  given  y 


y^l\  we  can  calculate 


. In  oilier  words. 


( ) (,,) 

given  as  an  initial  val.ua  we  can  calculate  ; , ^hen,  using  tr 

recurrence  relation  (3-8)  we  on  calculate  y from  wiiich  we  now  can 
calcul  ite  . Tiiis  recurrence  procedure  can  then  be  continued  for  a., 
many  I n vs  as  desired. 

A power  series  solution  of  these  types  of  models  is  usually  "absolutely 
convergent"  in  the  complex  variables  sense.  "Numeric  convergence"  of  a 


power  series  with  finite  arithmetic  is  not  absolute.  For  tills  reason  the 


recurrence  relation  will  l>c  "numerically  good."  only  v/idiin  a given  r.ei/.h 


NUMERIC  CONTINUATION 


f 


1 

- 

i 


J 

l 


i 


« 

! 

> 

i 

I 
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The  power  series  representation  of  the  solution  of  the  scale!'  integro- 
differential-difference  equation: 

t 

y(t)  = f(y(t))  + / g(t-X)  y(x)  dX  (-1—1) 

n 

h..:>  1a.-.!:;  shown  to  be  amenable  to  calculation  with  a reasonable  effort. 

The  rower  series  has  a finite  range  or  neighborhood  about  the  center  of 
expansion,  within  which  tie  solution  values  (as  contrasted  with  the  to'.-.-r 
series  coefficients)  can  be  calculated  with  sufficient  accuracy’. 

V-.e  will  denote  by  T the  range  of  the  independent  variable  t e[o,i'3 
which  yields  y(t)  with  sufficient  accuracy.  Ke  will  call  the  process  of' 
obtaining  new  coefficients,  about  t = T,  numeric  continuation . It  is  quite 
siiailar  to  analytic  continuation  from  complex  variable  theory.  In  tins 
process,  we  use  the  acceptably  accurate  values  y('i)  as  initial  values  and 
generate  the  additional  power  series  coefficients  by  use  of  recurrence 
relations.  This  is  well  known  in  ordinary  differential  equations,  but,  a 
slight  problem  must  be  handled  for  the  .integro-difference  term.  This 
problem  arises  from  the  rianory  or  history  contributions  of  the  intogro- 
dif fere: ice  term. 

introduce  :,  a shifted  independent  variable,  t = t - T.  Ke  re- 
writ  -•  ■ ' 1)  as: 

T T+  t 

;/(T+r)  = _'(y(T+:)}  ■ j g(T+t-X)  y(\)  ■ J g(?+i-\)  . i\. 

o T 

(4-2) 


9 


1 


t 

i 

1 


J 


,7 


(•1-10) 


11 


By  inspection  we  f Lnd  all  values  in  the  above  equation  arc;  known, 

{ ■ ) 

is  a sir, pile  process  to  calculate  each  c, 

We  can  also  represent  the  second  integral  term  in  (1-2),  as  a 
series: 

/ y(\)  dk  = J /:  \ 

T i=a 

shifting  the  1 units  of:  integration  giv. 

T+\  T 

j g(ihi-X)  I/O)  dX  = J u (!->.}  ; ,:\)  dX 

X O 


T:  .a  subscript  T indicates  that  those  power 


series  coefficients  ar 


the  center  of  expansion,  t=T . Coinparing  equation  (4-12)  with  (2- 
' but  the  'n,g  ^ is  calculated  like  (2-1J) 


(i) 


(J)  „ (i-1- 


«/ 


(i-Dl 

(i+j-1)!  ' 


a:  id 


Jo) 


0.  These  ccefficients  are  : r the  shifted  iivd.p-..  d nt 


The  power  series  in  terms  of  the  shifted  independent,  variabl 
;o  bo  numericcrlly  ••wild  over  a finite  ■.  Further  continu'd, 
.res  would  ch  u.  :•  • s'  nificantly  • he  expense  of  slot 

as  over  all  pr  .v.io’  ,s  ranges. 

Via  have  ps.  rented  . hods  to  caLcul.r  ■ ",  n.tctly": 
i T i 

j / d ' ’ - dk  + j 

o o o 


- I 


<1, 


(i) 


(i)  i 


-ore 


+ i . 


and  it 
power 

(4-11) 

(4-12) 

e about 
4)  we  see 

(4-13) 

variable 

e will 
ion  pro- 
pov.er 


l X)  dX 

(4-14) 


+ 


(4-15) 


12 


The  usual  form  of  Use  kernel  y(x)  is  the  negative  exponential,  y ?::)  -- 
■ :r(-x).  See  Nachlinger  and  V'Jheeler  (1971) , Nunziato  (1971),  Cushing 
(!  976).  Assuming  this  form,  we  rewrite  lie  first  integral  term  on  the 
right  side  of  (4-14): 

T T 

/ exo(-(t+>->:  > j\)  ja  - / -.-i  t:Xp(-( t-a;;  //(a; 

o o 


l 

- ecyC-.l  j - .(-(x-X))  ;.(X)  . (4-16) 

a 

i 

Vvo  note  that  the  contribution  of  ./  (iron  equations  (4-3)  through 

(4-1 0) ) is  like  a "damr-U"  . The  • contribution  s’.tould  be  . .rail 

p 

e : pored  to  Ji  ^ if  the  value  of  T is  large  ■ rough  to  cause  sig:.i;  Leant 
d a ping  exp(-T)  . Previous  experience  in  ord  i n. iry  differential  oyn  ions 
'•  iicate  that  large  ranges  are  probable,  using  say  twelve  to  twenty  terms  of 
a.'  tower  series  Do  iron  (1970),  Hunter  ga.i  Childs  (1977).  U’e  car.  shift  a 
power  series  to  an  approximate  scries  at  the  next  center: 


l 

i-o 


x 

r . 


co 


.1 


(i) 


(4-17) 


n+1 


(,’  , . - T ) = t - A , 
r.  <-x  n n > rl 


(4-13) 


n+J 


(■)  v-  (k)  / k \ . 

' .)■.  h”  (i}  \<! 
n 


the  proper  power  . ;e i -s  substitutions  e.:ut‘.  Lon  (4-2)  boo.  *.• 

0 - ci> 

1 <i+u  . . r;/;j  / - i u 


rj  ■ , v a)  ; 

i-o  ‘ i.  i~o  ' >i 


(4-19) 


(i)  i 

1 


(4-20) 


•••uich  yields  the  i -currence  relation: 


13 


(i+1)  ij 


(i+1) 


(i) 


CIT 


(-:■ 


+ r.r 


(i)  i 
1 


OL- 


(i+1) 


1 

i+1 


L (i) 
n 


+ 


(4-21) 


Cur  c’.;.-'itpitational  strategy  is: 

(1)  input  the  initial  value  and  initialise  the  q coefficients 

to  zero.  Note  that  tit. ■ first  center  of  is  . . 

(2)  t> j now  calculate  the  a coefficients  about  tie  current  center 
of  expansion.  If  the  current  center  of  ■.o-punsion  is  T , , then 
use  equation  (4-10)  to  calculate  the  .-  c uficicnts.  If  the 
current  center  of  expansion  is  greater  than  T^,  then  use  tin- 
shift  algorithm  ( 13)  and  (4-19)  to  calculate  the  q coof:  i- 
cients  about  the  center. 

(2 ) ( • ) 

(3)  Calculate  the  coefficients  y •••  .'  by  the  recurrence 

relation  (4-21) . 

(4)  Determine  tire  range  of  the  power  series  v.hich  will  yield 
specific  accuracy  over  the  interval  d.  .aired.  Then  wa  nay  new 
evaluate  the  ■ ( ) at  any  point  within  the  range  calculated. 

( ,)  if  the  rr.  e does  .. >t  extend  over  the  h larval  desired,  V'\\ 

evaluate  the  it.rcri.cn  y at  the  c.  ..  e :•  f.t  of  the  range  .on 
becomes  ti: : next  ..cor  of  exponsroa.  - at  the  cor.'ri  .-  c 
process  beg i nr.'  :j  v i cii  step  (2). 


CHAPTER  V 


TiiE  1'OUNDARY- VALUE  PROBLEM  IN  INTEGKO- 
DlrTERENTIAIi-DIE'FEREJ.CE  EQUATIONS 


A general  fora  of  .linear  integro-dif f e c- ntial-dif ference  nq’.-itiQns 
re -Id  be 


y = Ly  + + j G(t-X)  y(X)  dX 


V'-i; 


.... ,_-re  I and  G arc  k ncwn  matrix  functions  of  In  addition  j and  j an 
vector  functions  of  i , and  f is  a known  vector  function  of  L.  A second 


ier  scalar  linear  integro-differentiul— dri-erence  equation  la.vo 


x + u x t >■;  = sin(t)  + J exp  (t-> > x(X)  dX 

o 


(5-2) 


is  a well  defined  problem  if  x(o)  and  x(o)  are  kncwn.  We  cai  convert  this 
:o  the  form  (5-1)  by  defining 
(1)  _ 


y 


y 


(2) 


* *| 

" 

L = 

-5  -v. 

| 0(a)  = 

0 

0 v:__  (-a) 

f(t)  ~ 

0 

. • ln{  t) 

_ 

- 

(5 

Most  higher  or:  r problems  in  integro-cif ferential-dif fcrence  equa- 
ls will  likely  be  of  the  same  type  as  this  example.  That  is,  one  oqua- 
t jn  In  a set  of  ordinary  differential  coin,  cions  will  have  an  integio 
. - f f erence  term. 

Tne  solution  of  boundary  value  problems  in  integro-dif ferential- 
difidrence  ec rations  of  a general  nature  are  possible  with  the  use  of 
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the  codes  o'.  Uoiron  (1971)  and  Hunter  and  Childs  (1977).  The  codes  must 
be  supplemented  with  subprograms  h o ! cmentations  of  conations  (2-12)  or 
(2-13),  equations  (4-9)  or  (4-10),  e.rjation  (4-13),  ar.d  equation  ’(4-19) . 


LIST  Of  RliduRENCES 


Acton,  F.S.,  Nu.iorical  Methods  That  work,  Mew  York:  Harper  and  Kow 
Publishers,  1970. 

Alien,  R.G.D. , Mi Cioi.utical  Analysis  for  j'conanists.  New  York:  St. 

Martins  Press,  19 so. 

Baker,  C.T.II. , "M.-thods  for  Inteyro-Dif  f.crontial  Equations" , in 
Numerical  Solui.  ion  of  Inte<.jral  Equation  ;,  .-viited  by  Delve:.,  L.M.  and 
Walsh,  J. , Oxford:  oxrord  University  Press , 1974. 

Bellman,  Richard  and  Cooke,  Kenneth  1..,  r r forential-Difl\ •i.vikv  E tua- 
tions,  New  York:  Ac.djuic  Press,  1963. 

Bellman,  Richard  and  Cooks,  Kennetii  L. , Meern  Elementary  Pi  ''..re: Rial 
Equations , Menlo  Park : Addison-Wesley  1'  .'  isiiing  Co.,  1971. 

Bownds,  J.M.  and  Wood,  Bruce,  "Oil  Numerically  Solving  Non-Linear 
Vol terra  Integral  Equations  with  Fewer  Computations",  SIAM  aw  anal 
on  Numerical  Analysis,  Vol.  13,  No.  5,  pp.  705-718,  1976. 

Cushing,  J.M. , "Forced  Asymptotically  i\  r iodic  Solutions  of  Predator- 
Prey  Systems  with  or  without  Ilereitary  Erf  acts",  SIAM  Jour  a 1 on  i 
Applied  Mr then-n tics , Vol.  30,  No.  4,  ip.  665-674,  1976. 

Cushing,  J.M. , "Periodic  Solutions  of  Voiterra's  Population  Equation 
witii  Hereitary  Effects",  SIAM  Journal  c.:i  Aooliad  Mathematics,  Vol.  31, 
No.  2,  pp.  251-261,  1976. 

Doiron,  H.E. , "An  Indirect  Optimization  Method  with  Improved  Conver- 
gence Characteristics",  Ph.D.  Dissertation,  University  of  K as Lon, 
Houston , Texas,  May  1970. 

Doiron,  II.E. , "Numerical  Integration  Via  Paver  Series  Expansi  on", 

M.S.  Thesis,  University  of  Houston,  Houston,  Texas,  August  1967. 

Feller,  W. , "On  the  Integral  Equation  of  Renewal  Theory",  Aim:  d of 
Mathematical  Hiaristics,  Vol.  12,  pp.  243-267,  1941. 

Hunter,  C.D.  and  Childs,  S.B. , "A  State  fical  Study  of  i ;.il 

Analysis  Ap;>.1.io.i  to  the  Regression  oi  nfa  Order  Different! •.  ; d nations", 
Transactions  of  the  TWenty-Second  Conference  of  Amy  Mi th  •:  ; -ions, 

Report  No.  77-1,  t.S.  Amy  Researcli  Office,  Research  Trie  . i > “ark. 

North  Carol in  1977. 


Kroysziy,  Erwin , Vivanced  Engineering  Mi fiiernatics , New  Yo.h:  violin  Wiley 
* Riio,  Lie.  , 1.96/. 

Lopidus,  !,.  and  Sienf  eld , s.n.,  Numei  leal  Solution  of  Ord.i t ury  Differ- 
ential Ecn.iat.ii  ns,  New  York : Academic  Press,  Inc. , 1971. 


17 


15.  Malone,  D.W. , and  Conner,  J.J.,  "Finite  Elements  and  Dynamic  Visco- 
elasticity" , Journal  of  the  Engineering  Mechanics  Division,  ASCE,  pp. 
11-15-1158,  August  1971. 

16.  Miller,  J.C.P.,  "The  Numerical  Solution  of  Ordinary  Differential 
Equations" , in  Numerical  /dialysis , edited  by  Walsh,  J. , Washington: 
Tbaapsoa  Book  Company,  1967. 

17.  Millar,  R.K. , "On  Vol terra1 s Population  Equation" , SIAM  Journal  on 
App.l io  i Mathematics,  Vol.  .1  !,  No.  3,  pp.  446-453,  In 56. 

16.  N.iehliu  ?r,  H. R.  and  V.a  -ah  L. , "A  Uniqueness  'i  . • ■ am  for  Rigid  . 
Heat  Co:  due  tors  with  i oi.y",  : irber.lv  of  A; -pi  f d t uvatics, 

Vol.  3.1,  Ho.  3,  pp.  267-273,  1373. 

19.  Nunziato,  J.W.  , "On  Heat  Cxmducrion  .in  Materials  \.vH'>  Memory" , 
Quarterly  ol:  Applied  Marie;:  a tics,  Vol.  29,  pp.  187-Hm,  1971. 

20.  Richtnryo.r , R.D.  and  Morten,  K.W.  , Difference  Methods  for  Initial- 
Value  Problems , New  York:  Join  Wiley  & Sons,  Inc.,  1367. 


21.  Roberts,  S.M.  and  Shipman,  J.S.,  Two-Point  l-u„  1 e y Value  Problems : 
Shooting  Methods,  New  Yo  k:  American  Elsevier  Pubiisliing  Co.,  1972. 

22.  Schui  l /. , U.J. , "Solving  Control  Problems  with  Discrete  controls" , 
a\i«Xk'r-ITC-02-08-76-019 , l.i-h.  635  Report  for  ir.tani  Training  Center 
in  Conjunction  with  Texas  AVM  University,  April  1976. 

23.  Scott,  Melvin  R.  and  Vn-.:  Vendor,  Walter  II.,  "A  Comparison  of  Several 

Invariant  Imbedding  Algor : . has  for  the  Solution  of  'IWo— Point  Bouraery 
Value  Problems",  Applied  v ueheiratics  and  Couputntion,  Vol.  1,  pp.  18/- 
213,  1975.  ~ ‘ """ 

24.  Scott,  Melvin  R.  and  Watts,  Herman  R. , "A  Systemic. d Collection  of 
Cod  -s  for  Solving  Two-I’osP.  Boundary  Value  ProbU  us",  SAND  75-0539, 
SanJia  laboratories , Alb’.*piurque,  New  Mexico  8/115. 


APPENDIX 


POKER  SERIES  ALGORI'i  PIS 
A Suimury 

Power  series  operations  are  straightforward  with,  the  use  of  auxiliary 
series  and  repetitive  applications  of  known  algorithms  for  series  addition, 
sui  .motion,  multiplication,  division,  and  elcsenfcary  functions-  Eh' 
follow?  eg  power  series  algorithms  (Doiron,  1970)  are  important  to  thin 


Addition:  f(t)  = g(  b)  + h(  t) 

f(i)  „ U)  „ h(i) 


plication:  f(t)  - f(j( ’),  h(t))  - g( t)*k( t) 

Jk)  v M , (k~i) 

J ~ L i i ,l 


ilvision:  f( t)  = j'(g(t)»  h(t))  = g(t)/k(t) 

Jo)  (o)  (o) 

1 = :i  /n 


~ I h" 

i-o 

Jo) 


(i)  (k~i) 

q 


k > 1 


Sine  and  cosine 


h(t) 

= CO,: 

(ll(D) 

1 

| 

g(t) 

= 6 in 

('j(t)) 

j 

■ (o) 

t: 

= COi'i 

(/°>) 

i 

(A-5) 

Jo) 

J 

= air, 

J(o)) 

(A— 6) 

Jk+I) 


X <■  y(i)  g(k-iJ 


/•  I v a 

v=o 
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